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Abstract 

A Bose-Einstein condensate of atoms, trapped in an axially symmetric har- 
monic potential, is considered. By averaging the spatial density along the 
symmetry direction over a length that preserves the aspect ratio, the system 
may be mapped on to a zero temperature noninteracting Fermi-like gas. The 
"mock fermions" have a state occupancy factor (>> 1) proportional to the 
ratio of the coherance length to the "hard-core" radius of the atom. The map- 
ping reproduces the ground state properties of the condensate, and is used 
to estimate the vortex excitation energy analytically. The "mock-fermion" 
description predicts some novel collective excitation in the condensed phase. 
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Recently there has been a renewed interest in the Bose-Einstein condensation(BEC) of a 
gas after its experimental demonstration with rubidium vapour in a trap at a temperature 
of 170 nanokelvin and at a number density p = 2.5 x 10 12 atoms per cc. This experiment 
has been followed by others using alkali atoms |2[. At these temperatures the atoms form a 
weakly interacting metastable gas of Bosons. For a non-technical account see the review by 
Burnett f| . In a typical device, atoms are trapped in a potential which is well described by an 
axially symmetric parabolic confinement. The oscillator frequency in the symmetry direction 
is larger than the frequency in the plane perpendicular to it. The experimental situation of 
interest to us is the one with rubidium vapour, where the s-wave scattering length between 
two atoms is known to be positive. The effect of the interatomic interaction may be mocked 
up by a repulsive pseudo-potential (4j. The interaction energy is propotional to ap 2 , where 
a is the s-wave scattering length and p is the number density of the atoms. The properties 
of the condensate have been studied by constructing the density functional involving this 
replusive interaction energy, and the potential energy of the atoms in the trap 0-0] • In this 
paper, we first note that by averaging the spatial density of the condensed bosons along 
one direction, it may be reduced to the same form as the density of a non-interacting Fermi 
gas. We chose the averaging direction to be the symmetry axis (the z-direction ), along 
which the harmonic confinement is steeper. This enables us to use the Fermi gas model to 
compute the low-lying planar excitations of the condensate in the shallow well. Moreover, 
the averaging distance is chosen to preserve the aspect ratio (the ratio of the length scales 
in the planar to perpendicular directions ) of the original trap. The Bose-condensate is now 
described by a three-dimensional noninteracting "Fermi" gas, trapped in the same planar 
parabolic potential as the original system, but free to move along the z-direction within 
the averaging distance. One peculiarity of these mock-fermions, as we call them, is their 
occupancy factor per state. Instead of being one (or zero ) at T = 0, it is multiplied by a 
large factor proportional to A^/2a, where = h/pp, and 2a the apparent size of the atom. 
In fact, the Fermi momentum hkp is such that kp 1 = £, where £ = (87rpa) -1 / 2 is just the 
coherance length in the bose condensate. It is remarkable that the kinetic energy of these 
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mock fermions exactly reproduces the condensate energy in the large-N limit. The latter 
is in fact calculated by neglecting the kinetic energy of the bosons. After having shown 
this equivalence, we go on to use this model to calculate some other properties. These 
include the velocity of sound, and the vortex excitation energy. The sound velocity in the 
mock-fermion ideal gas is the same as in the Bose-condensate. A simple estimate of the 
vortex excitation energy is made by "digging" a hole in the central density, i.e., promoting 
all the s-state mock-fermions out of the Fermi sea to states of non-zero angular momentum. 
This reproduces the numerical results of Dalfovo and Stringari || satisfactorily. The latter 
calculation involved solving a nonlinear Schrodinger equation that was obtained from the 
density-functional formalism. Finally, our description also predicts some novel collective 
excited states with zero angular momentum involving a large number of mock-fermions. 

We begin with the ground state energy for condensed bosons given by the Ginzberg- 
Gross-Pitaevskii || energy functional, 
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where m is the mass of the atom, u±,u>3 denote the oscillator frequencies in the transverse 
direction and in the direction of the symmetry axis( z-axis), and a is the s-wave scattering 
length which defines the strength of the interaction in the pseudo-potential method. The 
condensate wave function is usually denoted as 

ip(r) = \J p{r) expzs(r). (2) 

We work in the limit of strong repulsive interaction where the kinetic energy term can be 
neglected in the condensate phase. In the large N-limit we can then obtain the density p(r) 
by minimising (E — pN), where p is the chemical potential. This is the Thomas- Fermi 
expression for the density within the classical turning points : 

P(r) = 7^-?("iri + ^)]- (3) 
Ann a 2 

The density is zero outside the turning points. Note that when the scattering length a is 
positive, the chemical potential p is necessarily positive since p(r) > 0. This is the situation 



with the BEC of rubidium atoms. It is now straight-forward to calculate the particle number 
and the energy using the density given above. The particle number is given by, 

N =I^ = mr x &- (4 » 

where 



a± = J—; a 3 = J—, (5) 
V mwi v 171^3 

The integration limits are set by the turning points in the z- and the redirections. These 
may be obtained by first performing the z- integration for a fixed r±, and then allowing the 
latter to vary within the prescribed limits. The condensate energy is obtained by substituting 
the expression @ for the density in Eq.(|l|), and is given by 

As noted before, the kinetic energy term in Eq. ([]]) is neglected. These are known results. 
Before we go further, we note that the ratio E/N = 5/i/7. This is indeed the energy 
per particle in a non-interacting fermionic system whose single-particle density of states 
g(E) oc E 3 / 2 , obtained by convolving the state density in a planar parabolic potential, with 
the density in the transverse direction. 

We now take the crucial step by averaging the density in the direction of the symmetry 
axis over the a scale L 3 : 

Pa(r±) = J^j dz P ^ ' ( 7 ) 

We refer to p a as the mock-fermion density for reasons that will be clear soon. The length 
scale L 3 (— L 3 /2 < z < L s /2) will be fixed later by fitting the aspect ratio. A little algebra 
shows that the result may be written as 

pa{r±) = -Ptf (r±), (8) 
a 

where 
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and the dimensionless parameter a is given by 

L 3 a 



a: 



ita 
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The Thomas- Fermi expression for the density, Ptf, given by Eq. (||), is in fact the density 
of spinless fermions confined in a harmonic potential in the plane and free to move in the 
z-direction. Note from Eqs. (§, 0) that for as — > 0, the system becomes two-dimensional, 
and there is no condensation since p a — > 0. Similarly, the condensation density is depleted 
as the scattering length a is increased. Note that 1/a plays the role of the occupancy factor 
which is not necessarily unity. It is in this sense we call these particles mock-fermions. 

To proceed further we first compute the number of particles N and the condensate 
energy E, using p a {r±) as given by Eq. (§). Once again the number of particles is defined 
by 



N = L 3 d 2 r ±Pa (r ± ) 



L3 / 2/x ^ 5/2 



(11) 
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Substituting for a from Eq . (|10"D reproduces the particle number given by Eq.(f|). The total 
energy in fermion-like picture is given by, 



E[ Pa ] = - [d 3 r 
a J 



— r(f) + U(f)p TF {f) 



(12) 



where r(r) is the kinetic energy density of the fermions in the Thomas- Fermi approximation, 
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and U(r) is the confining oscillator potential in the plane. The particles are of course free 
to move in the z-direction, but confined within the length L 3 . Substituting for the density 
and the kinetic energy density in Eq . fll2T) , we immediately obtain the total energy 

L 3 . 2/i 



E 
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Again, substituting for a yields the energy of the Bose condensate as given by Eq.(§). This 
is remarkable since the energy of the interacting bosons in the BEC is identical to that of 
the non-interacting mock-fermions. 



Note, from Eqs. ( |TT| , that both the particle number N and energy E are independent 
of the choice of L 3 . This is because it comes in the combination L 3 /a. This does not mean 
that the choice of L3 may be arbitrary, since we shall presently see that that the excited state 
properties depend crucially on it. A physically meaningful way of fixing the length is to set 



it through the aspect ratio y < x 2 > / < z 2 > = lo 3 /uj±. Consequently the size of the cloud 
is approximately the same in BEC and in the mock-fermion picture. We have cylindrical 
symmetry in the latter instead of ellipsoidal symmetry of the BEC. In the fermion picture 
< z 2 >= L\jYl. The average distance in the plane is (< r\ >= is the same as in the 

BEC. Equating the aspect ratios in both pictures we get, 



24/i 



y 7mw 3 

This then determines the occupancy factor a uniquely, 



^ = \h=7S • (15) 



where a is the scattering length. Note that a now depends on the chemical potential. To 
obtain the sound velocity in the interior of the cloud, we set the density p a from Eq. (§) to 
its central value: 

1 f 2m I 3 / 2 

P« = T^r tt p ■ ( 17 ) 

But from Eq. (|16j), a is proportional to y/JI, yielding the correct linear depence between the 
density p and p as in BEC. We then obtain the same sound velocity u s as in BEC : 



2 p a Op p 



m dp a m 



(18) 



By contrast, u 2 = 2p/3m in a free Fermi gas with constant a. By setting k 2 F = ^pr, we see 
from Eq. (|Tj]) that a is proportional to the dimensionless quantity akp. This results in the 
occupancy factor 1/a being 
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Since the de Broglie wave length Xp = 2it/kp may be regarded as the resolving power, we 
may interpret the degeneracy described by the parameter 1/a as the collective number of 
the atoms which can be accommodated within a wave length. The quantum mechanical 
wave functions of all these atoms overlap substantially over a wave length. Thus a single 
"fermion" in our picture is as if made of 1/a number of mock-fermions. It is also easy to 
check that 



where £ is the coherence length of the bosons in the condensate ||. Using the parameters of 
the experiment [|l[] , it is straight-forward to get a numerical estimate of this collectivity. For 
example, for iV = 5000, we find 1/a = 140. Note that in our picture the mock-fermions are 
essentially free apart from the confinement in the plane and over a length in the z-direction. 
The single fermion excitations in this picture may now be described as the excitation of 1/a 
mock-fermions. This has interesting consequences in computing the excitation energies of 
the system. 

As an immediate application we may consider the vortex states discussed by Dalfalvo 
and Stringari || . In our Fermi gas picture, we may make an estimate for a vortex excitation 
in the static "shell model" picture. For making this estimate, consider the ground state in 
which mock-fermions multiply occupy the states of the two-dimensional harmonic oscillator 
upto the Fermi energy. Since the density profile of a vortex has a hole in the centre, we may 
simulate it by promoting all the s-state mock-fermions out of the Fermi sea to states with 
/ 7^ 0. Each of these s-states contain 1/a mock-fermions in our picture. Keeping N fixed, 
the number of / = states below the chemical potential fi is given by fi/(2huj) for large 
fi. Since the degeneracy of the harmonic oscillator state just above the the Fermi energy is 
(yU + l)/huij_, it is possible to promote all the s-state particles to this state (upto a maximum 
angular momentum). The excitation energy per particle in the presence of a vortex is then 




(20) 
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given by, 



AE/N 




2Tiuj_ L 



(21) 



This should be compared with the energy hQ c calculated by Dalfovo and Stringari |J for 



n — 1, where Vt c is the critical angular velocity and Tin the angular momentum of the vortex. 
In Fig.l, we display the result of our estimate for Q c , obtained from Eq. (^l|), with the graph 
given in || for the same. The agreement for large N is good between our calculations and 
those obtained by solving the non-linear Schrodinger equation. For comparison we have also 
displayed the calculated values of fl c obtained in the large N limit by Baym and Pethick || . 
While the power law dependence on the chemical potential is the same in our calculation 
and in ||, we cannot get the dependence on the logarithm of the chemical potential in the 
static shell model type of analysis. In our naive model, the vortex energy is insensitive to 
the vortex angular momentum Tin upto a maximum. 

As pointed out in the introduction, the mock-fermion picture also predicts collective 
particle-hole excitations from the filled Fermi level to an excited state, where all - mock 
fermions of a given angular momentum state are excited together. Since the adjacent shells 
of the harmonic oscillator have states of opposite parity, excitations to the next higher state 
of the - mock fermions whould entail a huge change in angular momentum. Such high-spin 
collective excitaions are unlikely. There could be, however, collective excitations of energy 
-huix and zero angular momentum. Both the above types of excitations are of a novel 
collective type, and involve a hundred or more mock-fermions. As such, they are much 
larger energy excited states than the ones predicted and analysed by Stringari || recently. 

To conclude, the Bose-Einstein condensate of interacting bosons is described in 
this paper by a noninteracting Fermi gas. We call them mock fermions since they obey a 
generalised Pauli principle with the occupancy of a state i given by 



1 

Tij , £j <~, fl , 



a 



0, €i > /i . 



(22) 
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This occupancy factor is directly proportional to the coherance length £ = (87rap) l / 2 divided 
by the "size" 2a of the atom. The collectivity is generated by the combined motion oil/ a 



mock fermions. In passing, it is worth mentioning that an occupancy factor of the type ([22 ) 
at T = arises naturally in systems that obey the so-called Haldane statistics [[UJ] . It has also 
been noted that in such systems that the single particle density is scaled by the occupancy 



factor ||TT|. Irrespective of this aspect, the present analysis of the bosonic condensate in 
terms of mock-fermions is intersting in its own right, although one should only expect it to 
hold for global properties, and not for the analysis of correlations. 

We gratefully acknowledge comments from Diptiman Sen and discussions with Jim 
Waddington. One of us (MVNM) thanks the Department of Physics and Astronomy, Mc- 
Master University for hospitality. This research was supported by National Science and 
Engineering Council of Canada. 



9 



REFERENCES 

[1] M.H. Anderson, J.R. Ensher, M.R. Mathews, C.E. Wieman, and E.A. Cornell, Science, 
269, 198 (1995). 

[2] C.C Bradley, C.A. Sackett, J.J. Tollett, and R.G. Hulet, Phys. Rev. Lett. 75, 1687 
(1995); 

K.B. Davies, M.O. Mewes, M.R. Andrews, N.J. van Druten, D.D. Durfee, D.M. Kum, 
and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

[3] K. Burnett, Contemporary Physics, 37, 1 (1996). 

[4] K. Huang and C.N. Yang, Phys. Rev. 105, 767 (1957); 

T.D. Lee, K. Huang and C.N. Yang, Phys. Rev. 106, 1135 (1957). 

[5] G. Baym and C. Pethick, Phys. Rev. Lett. 76, 6(1996). 

[6] F. Dalfovo and S. Stringari, Phys. Rev. A53, 2477(1996). 

[7] A. Fetter, Preprint Stanford University, 1995 (unpublished). 

[8] V.L. Ginzburg and L.P Pitaeviskii, Sov. Phys. JETP 7, 858 (1958); 
L.P Pitaevskii, Sov. Phys. JETP 13, 451 (1961); 

E.P Gross, Nuovo Cimento 20, 454 (1961); J.Math. Phys. 4, 195 (1963). 
[9] S. Stringari, University Trento Preprint, March 1996(unpublished) 

[10] F. D. M. Haldane, Phys. Rev. Lett. 66, 1529 (1991); 67, 937 (1991). 

[11] Diptiman Sen and R.K. Bhaduri, Phys. Rev. Lett. 74 3912 (1995). 

R.K. Bhaduri, M.V.N. Murthy and M.K. Srivastava, Phys. Rev. Lett. 76, 165 (1996). 



10 



FIGURES 

FIG. 1. The critical frequency f2 c (in units of u±) is shown as a function of the particle number 
in the fermion-like picture (solid line). The dotted line is based on the graph given by Dalfovo and 
Stringari[6]. The dashed line is the large N prediction [5]. 
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